Quantum computing with a single molecular ensemble and a Cooper pair box 
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We propose to encode quantum information in rotational excitations in a molecular ensemble. 
Using a stripline cavity field for quantum state transfer between the molecular ensemble and a 
Cooper pair box two-level system, our proposal offers a linear scaling of the number of qubits in our 
register with the number of rotationally excited states available in the molecules. 

PACS numbers: 03.67.Lx, 85.25.Cp, 33.90.-|-h 
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One obstacle which transcends ah implementations of 
a quantum computer [TJ [2l |3l [4] concerns the extension 
of current proof-of-principle operations beyond a hand- 
ful of qubits. At the heart of this obstacle lies the expo- 
nential scaling of the Hilbert space dimension with the 
number of qubits. If one chooses to work with qubits 
encoded in separate particles, the available state space 
is exponentially large in the particle number, but se- 
lective access to individual particles and precise control 
of the interactions among individual quantum particles 
presents a formidable challenge. This has spurred inter- 
est in quantum systems that intrinsically support a vast 
Hilbert space. Obvious candidates are molecular quan- 
tum systems which easily provide 100 accessible internal 
rotational and vibrational levels O [6]. The quantum 
information capacity of such systems corresponds, how- 
ever, to a mere log2(100) ^ 6 qubits, and most molecular 
implementations to date have not exploited the rich inter- 
nal structure, but have focussed on other advantages pro- 
vided by molecular systems such as the large intermolecu- 
lar dipole-dipole coupling [7 , switchable interactions [8^, 
and long coherence times [6 . These advantages also make 
molecules very attractive for hybrid quantum comput- 
ing schemes involving solid state, optical and molecular 
quantum degrees of freedom simultaneously. Notably, in 
[9], it has been proposed to trap a mesoscopic molecu- 
lar ensemble at an antinode of the quantized field of a 
stripline cavity with a Cooper pair box (CPB) placed 
at the adjacent antinode. This setup is illustrated in 
Fig. [l^. The energy scale for the stripline cavity mode 
matches typical energies for rotational excitations of po- 
lar molecules, providing a natural interface between the 
cavity and molecular degrees of freedom. The large elec- 
tric dipole moment of polar molecules makes the strong 
coupling regime relatively easy to achieve while strong 
coupling of the field to the CPB has been demonstrated 
experimentally in [10||TT]. Furthermore, by using an en- 
semble of N molecules one achieves a ^/N enhancement 
of the coupling to the weak quantum field compared to 
the single molecule vacuum Rabi frequency g. In [9 , the 
essential idea is to counteract the rapid decoherence in a 
Cooper-pair box by transferring the quantum state to the 
molecular ensemble for storage of the qubit in a collective 
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FIG. 1: (Color online), (a) A Cooper pair box is strongly 
coupled to a stripline cavity as realised in ^OJ. Additionally 
a cloud of cold polar molecules interacts with the quantized 
cavity field, (b) Internal level scheme for a single molecule. 
The reservoir state |0) is selectively coupled to each of the 
excited states \i) by Raman transitions, (c) Encoding of K 
qubits in the symmetric states of an ensemble oi N > K 
identical particles. 



molecular excitation between quantum gates. 

In this paper we shall present a method for many-qubit 
quantum computing with a single molecular ensemble 
and a Cooper pair box. We shall apply an ensemble 
of N cold polar molecules with a ground state and K 
accessible excited states as illustrated in Fig. [TJd. The 
potentially available Hilbert space for the molecular sys- 
tem is of dimension {K + 1)^, but, by limiting ourselves 
to the symmetric states with at most one molecule pop- 
ulating each of the excited states, this is reduced to 2^, 
the Hilbert space dimension of a i^-qubit register. 

We assume all molecules are initially prepared in the 
ground state |0), which is coupled to the excited states 
\i) through a Raman process involving the cavity field 
coupling constant g and a classical field Vti(t). Since 
both fields couple symmetrically to all molecules in the 
ensemble, elementary excitations produce the symmet- 
ric (Dicke) states \i) = (I/VN) Y^. IO1O2 ... z^- ... Oat) 



and so forth, where the index j runs over ah molecules 
in the cloud. We define the collective raising operator 



(l/vTV) V . \i)jj{0\. In the regime with only few 
excited molecules, the collective operators approximately 
obey the bosonic commutator relation [mi^ml] ~ 1 and 
the cloud can be treated as a collection of K uncoupled 
harmonic oscillators. 

The conventional approach to encoding qubits in atoms 
and molecules is to encode a single qubit in a single par- 
ticle [71 [12], or in a single collective degree of freedom 
in an atomic or molecular cloud [H [13]. However, as 
was recently proposed in [14 , one can encode K qubits 
in a single cloud of identical particles, each with K -\- 1 
accessible internal levels by associating the logical regis- 
ter state \a1a2 ■ - ■ clk) {cii = 0, 1) with the collective state 
Ili{ml)^^ IO1O2 . . . Oat) with a^ particles populating the i^^ 
level, illustrated in Fig. [ijc). The permutation symme- 
try among particles is important, and, for instance, the 
logical two-qubit state |01) is an entangled state with no 
particles populating the excited state |^ = 1) and a unit 
population of state \i = 2) evenly distributed over all the 
molecules in the ensemble, in contrast to the conventional 
simple product state encoding of the same state. The ad- 
vantage of this encoding is that it circumvents the need 
for addressing of individual particles, since qubit access is 
granted by selective coupling of the ground state to one of 
the K excited states. The complication of the method lies 
in the restriction of the dynamics to the specified state 
space with at most one particle populating each of the 
excited states, and in the operations on these states that 
depend on the population of the other excited states. In 
[14^ it was proposed to use the Rydberg blockade mech- 
anism for controlled dynamics of neutral atomic ensem- 
bles. In this Rapid Communication we shall describe how 
the cavity field and the two-level system offered by the 
Cooper pair box can be used to achieve the same goal for 
polar molecules. 

We now turn to the setup indicated in Fig. [ija). The 
Cooper pair box is a superconducting circuit with an is- 
land onto which charge may tunnel through an insulating 
barrier as described by a phenomenogical Hamiltonian, 
H = -Ej/2^^\n){n^l\ + |n + l)(n| [15]. Due to the 
quadratic nature of the electrostatic interaction the en- 
ergy levels are non-equidistant, and using resonant tran- 
sitions only, the system may at cryogenic temperatures 
be restricted to the two lowest quantum states with cor- 
responding raising and lowering operators, a~^ and a~ . 
The CPB is mounted on a superconducting stripline cav- 
ity which can hold a cavity field with creation and annihi- 
lation operators c^ and c with very modest field damping 
[16]. The combined system CPB-cavity system is gov- 
erned by the Jaynes-Cummings type Hamiltonian 

HcPB = gc{cr~C^ + cr+c) + dcPB{t)cr^cr~, (1) 

where (^cpb(^) = ^cpb(^) — ^c is the tunable CPB de- 
tuning with respect to the cavity field. A number of phe- 



nomena related to the Jaynes-Cummings Hamiltonian in 
quantum optics have been observed in the CPB-cavity 
system [TOl [11] , and two-qubit gates on two Cooper pair 
boxes coupled to a single cavity field mode have recently 
been demonstrated ^. The molecular ensemble is ad- 
dressed by a Raman transition involving the cavity field 
and a classical field with tunable frequency and real am- 
plitude Qi{t) [see Fig.jljb)]. In the rotating wave approx- 
imation after adiabatic elimination of the excited state 
|e) the coupling of the cavity field to the z'th molecular 
qubit is described by 
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gi(t)(miC^ + m\c) + 5i(t)m\mi. 



(2) 



Here gi{t) = Qi{t)g\^7%/2A is the effective coupling 
strength with A the detuning with respect to the inter- 
mediate excited state, and Si{t) is the two-photon Raman 
detuning of level |i), cf. Fig. [lib). Coupling by higher 
order Raman processes with the cavity field and multiple 
classical field components allows exploration of a wider 
range of molecular states for which Eq. Q applies with 
modified expressions for gi{t). 

We will now describe how quantum information can be 
encoded and processed in the combined system. Initially 
the molecular ensemble is prepared with all molecules 
in the zero state corresponding to all qubits set to the 
value 0. We now need to specify how to carry out re- 
liable one- and two-bit gates on the system, and specif- 
ically for the ensemble encoding, we have to ascertain 
that no register state is populated by more than a sin- 
gle molecule. The cavity and CPB are also prepared 
in their ground states, and SWAP operations of arbi- 
trary unknown states between any qubit component of 
the molecular memory and the CPB via the cavity field, 
combined with an arbitrary single qubit rotation of the 
CPB two-level system by resonant driving with a classical 
field, implements this rotation on the desired single qubit 
of the register. As a fully entangling two-qubit gate we 
propose to SWAP one molecular ensemble qubit to the 
CPB and then SWAP another molecular ensemble qubit 
to the cavity field. The CPB-cavity interaction can then 
provide a state dependent phase, remaining with the two- 
qubit states when they are finally returned to the collec- 
tive molecular ensemble states. Before presenting details 
of these processes we note that the main elements of these 
steps are, indeed, very similar to the ideas for quantum 
computing with atoms coupled via a cavity field [17] and 
with trapped ions, coupled to each other via their collec- 
tive motional degree of freedom [18] • But we emphasize 
the significant difference that our molecules do not need 
to be individually addressed, and that we need to pass 
the cavity excitation through the Cooper pair box to re- 
strict the Hilbert space to two states per qubit degree of 
freedom and to provide the interaction in the system. 

The SWAP operations can be realized by adiabatically 
sweeping the detunings across resonance. When transfer- 
ring a molecular state to the empty cavity the coupling 



gi{t) is turned on with 6i{t = 0)/^^ ;» 1. As 6i{t) passes 
through resonance each basis state adiabatically evolves 
into the corresponding dressed state 

|+,n) = cos(^)|l)^|n), + sin(^)|0)^|n + l), (3) 
|-,n) = -sin(^)|l)^|n)e + cos(^)|0)^|n + l)e (4) 



with tan(2^) = 2gi^/n + l/5i{t) and energies E±^n = 
\5i{t) ± \\/5l{t) +4^|(n + l). Here n + 1 is the number 
of elementary excitations i.e. n = for the states |0)^|l)c 
and |l)m|0)c, and n = 1 for |l)^|l)c. As a result each 
state acquires a nonlinear phase relative to the ground 
state of (^±,n = \Jq -^i(t) T ^5^,(t) ^ 4.g^,{n ^l)dt. 
Thus the logical state |1) acquires a phase relative to |0) 
while following |+,0) (see Fig. |2J left panel). When re- 
turning the state to the ensemble the |l)-component can 
either be made to follow |+,0) (by sweeping from nega- 
tive to positive detuning) or | — , 0) (by sweeping from pos- 
itive to negative detuning). By choosing the latter option 
the total dynamical phase exp{— i J^ [^+(^) + E-(t)]dt} 
is exactly cancelled for any sweep that is antisymmet- 
ric about t = T/2. The |1) state then only acquires 
a geometric phase of e~*^. The effective coupling of the 
ground to any particular excited states is amplified by the 
number of participating molecules, i.e., it is g\fN^ with 
A^o the number of ground state molecules. This number 
has a quantum mechanical uncertainty when some of the 
other excited states are in superposition states of being 
populated and unpopulated. The relative variation in 
the coupling strength due to this uncertainty can be sup- 
pressed by requiring N ^ K^ but we observe that since 
we accumulate opposite phases in the two adiabatically 
swept passages, indicated by the dotted {N > Nq) and 
dashed {N < Nq) curves in the figure, the difference will 
effectively cancel. 

The phase dynamics is most easily visualized by view- 
ing the dynamics as a spin- 1/2 precessing about a ficti- 
tious magnetic field B{t) = —(^^,0, ^6i{t)). Two consec- 
utive sweeps of Si{t) in the same direction corresponds 
to a 27r rotation of B{t) giving a geometric phase of 
g-zTT ^]^ig SWAP operation thus incurs an extra single 
qubit Z-gate which must be absorbed into future oper- 
ations. When transferring between the cavity and CPB 
we have an always on interaction gc and two consecutive 
sweeps must be in opposite directions (see Fig. ^ right 
panel). This corresponds to following the |+,0) (| — ,0)) 
state twice and hence there is no cancellation of the dy- 
namical phase. However in this case the geometric phase 
vanishes since the fictitious B-field traces out a path en- 
closing a vanishing solid angle. The sweep may then be 
chosen such that the dynamical phase becomes mod 27r 
since there is no fiuctuation in the coupling ^c- For in- 
stance the cubic parametrization ^cpb(^) = So{2t/T — l)^ 
with 5o = 19.24^c and T = 20.77^"^ gives cancellation of 
the dynamical phase. The cubic form was chosen since 
it is odd about t = T/2 and offers faster implementation 
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FIG. 2: Left: In order to transfer a molecular state to the cav- 
ity we turn on a far-detuned field with bij gi ^ 1 and sweep 
through resonance adiabatically following the dressed state 
|+,0). To transfer the state back to the molecular ensemble 
the field is again turned on at Si/gi ^ 1, this time follow- 
ing the | — ,0) state through resonance. Phase errors due to 
fluctuating A^o (illustrated by the dashed and dotted lines) 
are exactly cancelled leaving only a geometric phase of e~*^. 
Right: When transferring states between the cavity and CPB 
one must follow the same path back and forth hence the ge- 
ometric phase vanishes. Since there is no fluctuation in the 
interaction strength gc the dynamical phase can be tailored 
to mod 27r. 



than a linear chirp while remaining in the adiabatic limit. 

To implement a single qubit gate on qubit i we transfer 
the molecular state to the cavity by tuning Si{t) across 
resonance and then transfer the state to the CPB by a 
similar sweep of ^cpb(^)- A single qubit gate may then 
be implemented on the CPB using microwave pulses [19] 
whereafter the state is transferred back to the ensemble 
via the cavity. To implement two-qubit entangling oper- 
ations we make use of the cavity-CPB coupling described 
by Eq. (fTl). The control qubit is transferred to the CPB 
and the target qubit is subsequently transferred to the 
cavity field as described above. From an initial detun- 
ing of ScpB/gc ^ 1 the CPB detuning is tuned close to 
resonance and back, the computational states evolving 
adiabatically along the dressed states as seen in Fig. [3] 
In order to implement a fully entangling controlled phase 
operation the functional form of Sqpb (t) must be chosen 
such that (p±^o = mod 27r and (p-^^i = n mod 27r. We 
find that parametrizing ScpB{t) = a(2t/T — 1)^+6 with 
a = 33.05^C5 b = 0.6664^c and T = 58.07/^c produces a 
controlled phase gate with near unit fidelity. 

Let us now address the feasibility of our proposal 
with current physical parameters. For the Cooper pair 
box, the dominant source of decoherence is second or- 
der charge noise, which is minimized by operation at the 
so-called sweet spot leading to dephasing times of the 
order T2 ~ 1 /is [20]. By comparison the vacuum Rabi 
frequency of the CPB-cavity field coupling is of the or- 
der gc ^ 27r X 50 MHz, so for single qubit rotations on 
the CPB we obtain ^c^2 ^ 300; i.e., one can implement 
on the order of 300 CPB-cavity SWAP operations be- 
fore the qubit decoheres. The conditional phase gate we 
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FIG. 3: A conditional phase shift is performed by exploit- 
ing that the dressed state energies are nonlinear in the ex- 
citation number n. Starting at fcpB(0)/^c ^ 1 the CPB is 
tuned close to resonance (fcpB(T/2) r-u g^ and back to end 
at fcpB(r)/^c ^ 1- The functional form of fcpB(t) is cho- 
sen such that (/?+,! = TT mod 27r while (/?±,o = mod 27r. 
The dashed line represents the dressed state | — , 1) which at 
^pb/^'c ^ 1 corresponds to |2)c|5')cpb i.e. a doubly excited 
state which is not part of the computational Hilbert space. 



have proposed has a gate time of T = 58.07^"^ giving 
T/T2 ~ 0.1. This ratio can be improved by more than 
an order of magnitude by replacing the conventional CPB 
with a recently improved so-called transmon design [20] , 
which operates in a regime where the dominant decoher- 
ence process is relaxation with T\ ~ 16 /is. The coher- 
ence time of the molecular ensemble is limited mainly by 
the collision rate and by coupling due to the long-range 
dipole-dipole interaction. In a magnetic trap, at T = 1 
mK the scattering rate due to the asymptotic r~^ inter- 
action is estimated to be 7 < 27r x 700 Hz [9 . In an 
electrostatic trap the induced dipole moment /iind leads 
to an T~^ interaction, but at T = 1 mK and /iind < 1 D, 
the above estimate for the scattering rate still holds and 
we conclude that the interactions within the molecular 
ensemble should not significantly lower the effective de- 
coherence time. With photon loss rates down to 27r x 10 
kHz and realistic values of gi up to 27r x 10 MHz [9], 
with current technology one could implement hundreds 
of gates, sufficient to provide proof of concept for the 
present scheme and to carry out simple error correction 
algorithms. 

In conclusion we have described a system which with 
current technology could provide quantum computation 
with an appreciable number of qubits. The ensemble 
encoding we have used is applicable to systems which 
contain a saturable element such as the Cooper pair 
box which furthermore serves to provide an entangling 
conditional phase gate. The molecular ensemble pro- 
vides an efficient encoding scheme for many easily ac- 
cessible qubits and long coherence times which may be 
improved even further by adopting a crystalline ensem- 



ble in which collisions are suppressed. If additionally one 
chooses a trapping scheme which is independent of the ro- 
tational state, there exist certain "magic" configurations 
for which detrimental coupling to the phonon spectrum 
is completely suppressed ^21^. Potentially in future ap- 
plications several ensembles could be coupled to an array 
of Cooper pair boxes in several interconnected cavities 
providing scalability as well as the benefit of parallel pro- 
cessing. Further improvements on the present work could 
include the application of Optimal Control Theory [6 to 
improve the gate time for the conditional phase shift as 
well for the SWAP operations. 

This work was supported by the European Commission 
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